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Abstract 

The rising cosmic ray positron fraction reported by the PAMELA collaboration has lead to a great 
^ deal of interest in astrophysical sources of energetic electrons and positrons, including pulsars. In 

^ this paper, we calculate the spectrum of synchrotron emission from electrons and positrons injected 

^-H from 376 young pulsars (< 10^ years) contained in the ATNF catalog, and compare our results to 

observations. We find that if objects such as the Vela and Crab pulsars have injected ~10^'^ erg or 
more in energetic electrons and/or positrons, they are expected to produce bright and distinctive 
fZ^ features in the synchrotron sky. Intriguingly, we predict hard synchrotron emission from these regions 

^ of the sky which is qualitatively similar to that observed by WMAP. 
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I. INTRODUCTION 



In polar cap [U |2] and outer gap [3j models, electrons are predicted to be accelerated in 
pulsar magnetospheres. The resulting electromagnetic cascades can include photons that are 
capable of producing highly relativistic electron-positron pairs. Through such mechanisms, it 
is plausible that pulsars generate a significant fraction of the GeV-TeV electrons and positrons 
present in the cosmic ray spectrum [3V(T0]. 

Recent observations from the PAMELA experiment and the Fermi Gamma Ray Space Tele- 
scope (FGST) have generated a great deal of interest in possible sources of high energy cosmic 
ray electrons and positrons, including pulsars |1H6]. In particular, the PAMELA collabora- 
tion has reported that the cosmic ray positron fraction (the ratio of positron to electrons- 
plus-positrons) increases rapidly between approximately 10 GeV and 100 GeV [TT] (see also 
Refs. [121 [13] )• This is in striking contrast with the predictions of models in which the cos- 
mic positron fiux is dominated by the secondary products of hadronic cosmic ray interac- 
tions [m Us]. More recently, the FGST collaboration has published their measurement of the 
cosmic ray electron (plus positron) spectrum between approximately 20 GeV and 1 TeV [16]. 
This measured spectrum contains a greater fiux of electrons above ~100 GeV than is predicted 
in typical galactic cosmic ray models [17], although not so much as previously reported by the 
ATIC collaboration [T8|. 

Among other plausible interpretations of these observations, it has been suggested that a 
nearby and relatively young pulsar (or pulsars) may be the source of the excess particles [H [3] . 
Assuming that the responsible astrophysical source is not unique in the Milky Way, other 
regions of our galaxy are also expected to contain large fiuxes of highly relativistic electrons 
and positrons. Such particles generate synchrotron photons through their interactions with 
the Galactic Magnetic Field, as well as inverse Compton emission by scattering with starlight 
and other radiation. Interestingly, emission consistent with synchrotron and inverse Compton 
emission from the inner Milky Way has been observed in the data of WMAP [T9| and FGST [20] , 
respectively, leading to the conclusion that a large fiux of highly relativistic electrons/positrons 
is present in the inner kiloparsecs of our galaxy. 

At this time, the origin (or origins) of the GeV- TeV electrons/positrons observed in the 
local and inner Milky Way has not been unambiguously determined. Among other sugges- 
tions, a great deal of attention has been given to the possibility that the excesses observed by 
PAMELA, FGST, and WMAP may be the products of dark matter particles annihilating [211 - 
|23] or decaying [21] in the halo of the Milky Way. In particular, if the dark matter consists of 
particles with a mass between several hundred GeV and several TeV, and annihilates or decays 
primarily to leptons, their products may be able to account for each of the aforementioned 
observations [211423] . 

In this paper, we consider the synchrotron emission from the electrons and positrons injected 
from nearby and young pulsars. If such emission could be observed and identified, it would 
provide a valuable test of the origin of the PAMELA and FGST excesses. In particular, the 
observation of synchrotron emission from regions around young pulsars (or the lack of such 
emission) would significantly weaken (or strengthen) the case for annihilating or decaying dark 
matter being the source of the PAMELA and/or FGST signals. 
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II. SYNCHROTRON FROM REGIONS AROUND YOUNG PULSARS 



Energetic electrons and positrons generated by pulsars move through the Galactic Mag- 
netic Field, steadily losing energy through synchrotron and inverse Compton scattering. These 
processes can be modeled by the propagation equation 



where due/dEe is the number density of electrons and positrons per unit energy, hE^ is the en- 
ergy loss rate, DE'^ is the energy dependent diffusion coefficient, and Q{Ee, x) is the source term 
which describes the spectrum and time dependence of injected electrons/positrons. Through- 
out this study, we use the following values for the propagation parameters: a = 0.43 and 
DE^ = 2.92 ■ lO^^cm^ s~^(£^e/GeV)", which are consistent with the primary-to-secondary ratios 
observed in the cosmic ray spectrum [20] ■ In estimating the energy loss rate, we adopt a mag- 
netic field density of pB ~ 0.6 eV/cm^, and a density of starlight and dust emission given by 
Pstariight+dust ~ 0.7 cV/cm^. Combined with the cosmic microwave background, these densities 
lead to 6 = 1.56 ■ lO'^^ GeV-^ s-\ 

For the source term, Q{Ee,x), we assume the injection spectrum to take the form of a 
power-law with a spectral index of a = 1.5, and exponentially cut-off above 600 GeV: 

f{z) = z--e-\ (2) 

where z = E/Ec and Ec = 600 GeV. This spectral shape is motivated in part by a desire to 
accommodate the rising positron fraction observed by PAMELA [4] , and by the pulsar emission 
model of Ref . [9] . The source term for an individual pulsar is then given by 

QiEe,x) = ^fiz)git)5'ix), (3) 

where Etot is the total energy injected by the pulsar in electrons/positrons, and g(t) describes 
the time dependence of the injection (normalized such that J g{t)dt = 1). 

The energy budget for a pulsar is simply given by its initial angular momentum, 

where Ms and Rg are the mass and radius of the neutron star, and Qo is its initial period. 
It is generally expected, however, that only a relatively small fraction of a pulsar's total 
spindown power will be transferred into energetic, escaping electrons/positrons. If ~10% of 
the spindown energy budget is expended in this way by a young and nearby pulsar (such 
as B0656+14 or Geminga), it could potentially accommodate the rising positron fraction ob- 
served by PAMELA Observations of the Vela pulsar also suggest that this is a plausible 
estimate [27] . With these motivations in mind, we adopt Etot = 10^^ erg as our canonical 
estimate throughout this study. 

The function g{t) accounts for the time dependence of the injected electron/positron spec- 
trum. We model this using two different functions. In the first model, we assume that all of 
the electrons and positrons are injected at the birth of the pulsar. Alternatively, we consider 
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the case in which the electrons/positrons are injected over an extended period of time during 
which the pulsar transfers its angular momentum, or spindown energy, into cosmic rays. In 
this model, the luminosity in electrons/positrons scales with g(t) oc (1 + t/ro)~^, where tq is 
the characteristic spindown time, given by jl] 

3c^I . /2xl0i2Q\2/ ^ \/15km\V Pr 



where B is the magnetic field at the surface of the star. In our calculations, we adopt a 
representative value of tq = 10^ years. In the remainder of this paper, we will refer to the two 
models of g{t) as the instantaneous injection and spindown models. 

To calculate the spectrum of electrons and positrons as a function of time and the distance 
from a pulsar, we have used the Green's function as given in Ref. |2H]- For the case of in- 
stantaneous injection, we were able to obtain an analytic expression for the electron/positron 
spectrum: 

= 3 ^"») [^(^' ^rn)]^^"^ exp [-X'^h{z, Xn,)] , (6) 

where h{z, x) = 2;^""/ (1— a;"~^), and = 1/ (1 — zC,). We have adopted the scaling parameters 
^ = bEcT, and A = t/vq, where = ADE'^j [(1 — a)hE^, r is the age of the pulsar, and r is the 
distance to the pulsar. 

For the spindown model the integral that gives the electron spectrum is given by 

To calculate the spectrum of synchrotron emitted by these electrons/positrons, we need to 
calculate the electron column density over a line-of-sight, dae/dEe. From Eqns. |6] and [7| we see 
that in both models the distance dependence is found within the same A^ term inside of the 
exponential. These integrals can be evaluated yielding error functions. The column density is 
thus found by making the following replacement in Eqns. [6] and [7} 



exp (— A^/i) — > exp (— A^/i sin^ 5) 1 + erf ^AdV^cos(5^ , (8) 

where S is the angle observed, and the distance to the pulsar, d, is related to the scaled distance 
by Ad = d/ro. 

At a frequency, u, the synchrotron power emitted by the electrons along a line-of-sight (per 
area, per solid angle, per frequency) is given by 

dP 27^3 6m> f^^^ , .da^ f mju \ 



du 1287r2 Jo dEe \z'^E^ub 

where ub = 3eB / {Anmec) ~ 21 {B/5fiG) Hz (we have used B = 5/iG throughout this study), 
and the function u{x) is given by [29j 



uix] 



da dyK,/^iy). (10) 

•Jo J xj sincf 



In the following section, we will apply this formalism to calculate the spectrum of synchrotron 
emission from the known nearby and young pulsars contained in the ATNF database. 
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III. RESULTS 



In this section, we consider known pulsars and estimate the synchrotron emission from the 
electrons and positrons they generate. In particular, we consider the pulsars contained in the 
ATNF pulsar database [32] , focusing on those which are relatively young, and thus most likely 
to generate the brightest signals on the sky. In particular, we include in our calculations the 376 
pulsars in the ATNF catalog that are less than one million years old. Electrons and positrons 
from older pulsars will have lost a significant fraction of their energy, and will have been diluted 
over large volumes of space, making them unlikely to provide a very bright or distinctive signal. 

In Fig. [T] we show contour maps of the predicted synchrotron emission from the known 
catalog of pulsars, calculated as described in the previous section. To aid in the comparison 
with the WMAP data, we have convolved the synchrotron spectrum with the sensitivities of 
the relevant WMAP detectors. For a given frequency band of WMAP, we calculate the specific 
intensity: 

f , dP , , 

Iu= diJ—w{u) (11) 

where w{h') is the sensitivity of the detectors normalized to unity. The maps in Fig. [l] show the 
signal in WMAP's 23 GHz band for the cases of instantaneous injection (top) and spindown 
injection over a timescale of tq = 10^ years (bottom). Overall, the two maps are quite similar, 
although they vary significantly in the regions near very young pulsars (such as the Vela or 
Crab pulsars). 

In Fig. |2| we plot the spectral index of the synchrotron emission from pulsars in our two 
models. To obtain the spectral indices shown, we calculated the 23 and 33 GHz maps in the 
manner described above and, assuming a power-law like form [I^, oc z^~"), we then calculate 
the value of a from the ratio the two specific intensities. Note that the spectral index, 7 (such 
that dN/dE oc E~'^), is related to a according to 7 = 2a + 1. 

At this point, a few comments are in order. Firstly, some of the most distinctive features seen 
in our contour maps can be associated with specific young and nearby pulsars. In particular, 
the emission from the region of the Gum Nebula (which contains the Vela pulsar, at I = —96.4°, 
b = —2.8°) dominates a significant fraction of the sky. A small region around the Crab pulsar 
(/ = —175.4°, b = —5.8°) is also clearly present. 

To understand why emission from the Vela and Crab pulsars appears so prominently in the 
synchrotron sky, consider the angular extent of the synchrotron emission from an individual 
pulsar, which is roughly given by 

».^.30"x»V r Y\ (12) 

d \ d y V 10, 000 years; ' ^ ' 

where Ddif is the diffusion distance scale, r is the age of the pulsar, and d is the distance to the 
pulsar. Of those pulsars in the ATNF catalog within 3 kpc (approximately 10 times the distance 
to Vela), there are only 6 pulsars with an angular extent of synchrotron emission smaller than 
produced by Vela. Of these, all but the Crab pulsar are considerably more distant, and thus 
less bright in synchrotron. Other nearby pulsars are bright, but too extended to provide a 
distinctive synchrotron signal. The Geminga pulsar, for example, is roughly as close to the 
Solar System as Vela, but is roughly 30 times older. As a result, its electrons fill a volume 
which encompasses the Solar System, making it difficult to detect its synchrotron emission on 
the sky, despite being one of the brightest sources of gamma-rays [31]. Furthermore, very young 
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FIG. 1: Maps in galactic coordinates of tlie synchrotron intensity in WMAP's 23 GHz band from 
the 376 pulsars in the ATNF database less than 10^ years old. The upper and lower frames show 
results assuming instantaneous injection, and a more gradual spindown injection (over a timescale 
of To = 10^ years), respectively. The plot uses a log-scale in Jy/sr with 0.1 between each contour, 
and with representative contours labeled (2.0 denoting lO^ *^ Jy/sr, for example). The normalization 
assumes that each pulsar generates 10^^ erg in electrons and/or positrons. See text for more details. 

pulsars other than Vela are in some cases too compact to distinguish from other emission from 
sources near the Galactic Plane. 

From Fig. |2| one can also notice that the spectral index of synchrotron emission from young, 
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FIG. 2: Maps in galactic coordinates of the synchrotron spectral index from the 376 pulsars in the 
ATNF database which are than 10^ years old. The upper and lower frames show results assuming 
instantaneous injection, and a more gradual spindown injection, respectively. In the upper frame, the 
contours run from 1.25 (dark regions) to 1.95 (light regions) in steps of 0.05. In the lower frame, they 
run from 1.5 to 1.9. 



bright pulsars varies with location. Around the Vela Pulsar (/ = —96.4°, b = —2.8°), for 
example, the spectrum is harder farther away from the source, which is a result of the fact that 
the most energetic electrons propagate away from the pulsar more quickly than its lower energy 
counterparts, leading to a particularly hard electron spectrum in the region of the surrounding 
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ring. 

Ideally, one would directly compare results such as those shown in Figs. [T] and |2] with obser- 
vations of the ~10-100 GHz sky. For example, we can compare our results to the observations 
of the WMAP experiment. Comparing our maps to the WMAP residuals (the emission not cor- 
related with known foregrounds) shown in Fig. 7 of Ref. [32], for example, we notice a number 
of features in common with our results. In particular, the most distinctive region of our maps 
(the region around Vela) is also one of the brightest in the WMAP residual maps (with resid- 
ual emission on the order of ~1 kJy/sr, which is not radically different from the ~0.2 kJy/sr 
found in this region of our map for a normalization of 10^*^ erg having been emitted from the 
Vela pulsar in electrons/positrons). Furthermore, this region possesses a notably hard spectral 
index, consistent with that predicted from cosmic electrons produced by the Vela pulsar over 
the past 11,000 years. In Fig. 7 of Ref. [32], one can also notice a small and bright region in the 
vicinity of the Crab pulsar (/ = —175.4°, b = —5.8°). As these represent particularly complex 
astrophysical regions, however, one should be careful to not draw overly strong conclusions 
regarding the nature of these residual emissions. 

The spindown time and the quantity of energy emitted into electrons and positrons is ex- 
pected to vary from pulsar-to-pulsar, potentially impacting the resulting intensity and spectral 
index maps. For those pulsars which we have found to provide the most distinctive signals (the 
Vela and Crab pulsars), however, our canonical values are appropriate estimates. Assuming a 
1.4 solar mass neutron star of radius 10-15 km, the 33 msec period of the Crab pulsar provides 
(2 — 5) X 10^^ erg of energy (and somewhat more if its initial period were shorter), with an 
estimated spindown time of several thousand years. In the case of the Vela pulsar, its observed 
period of 89 msec again leads to an estimated total energy on the order of 10^^ erg, and a 
spindown time on the order of 10^ years. 

The region in which our maps most significantly depart from those shown in Ref. [32] is 
in that surrounding the Galactic Center. The residual emission in this region, known as 
the "WMAP Haze" is likely hard synchrotron from a population of highly relativistic elec- 
trons/positrons. This is the topic of the following section. 

IV. COMMENTS ON PULSARS AS THE SOURCE OF THE WMAP HAZE 

The residual emission known as the 'WMAP Haze" has a number of intriguing charac- 
teristics. Firstly, it is consistent with being synchrotron emission from a population of elec- 
trons/positrons with very hard spectral index, dNe/dEf. oc E~"', 7 ~ 1.8. This is much harder 
than expected from supernova remnants (after the effects of propagation and energy losses are 
taken into account), for example. Efforts to identify the WMAP Haze with other possible emis- 
sion mechanisms (thermal dust, spinning dust, or soft synchrotron as traced by low frequency 
maps) have proven unsuccessful. Furthermore, as opposed to correlating to the disk or ISM of 
the Milky Way, the excess haze emission is approximately radially symmetric with respect to 
the Galactic Center, with an extent of ~20°. These characteristics have motivated the possi- 
bility that the haze emission is synchrotron from high energy electrons/positrons produced in 
dark matter annihilations taking place in the inner kiloparsecs of the Milky Way [22] . 

Alternatively, it has been previously suggested that a large population of (several thousand) 
pulsars in the inner galaxy may collectively produce the electrons and positrons leading to the 
observed synchrotron haze [33] . In that study, the authors argue that the very hard observed 
spectrum can be accommodated with an injected spectral index similar to that considered in 



8 



this paper, assisted in part by the diffusion hardening that occurs as a result of the propagation 
model. Furthermore, for appropriate choices of the model describing the spatial distribution of 
magnetic fields, they show that the morphology of the pulsars' synchrotron emission could be 
sufficiently similar to radially symmetric to avoid any obvious inconsistency with the observed 
haze [^). Note that our maps do not contain a particularly bright haze-like {ie. radially 
symmetric) signal, in part because our sample of known pulsars does not contain a large fraction 
of those located in the inner Milky Way. A more complete pulsar catalog would presumably 
yield a more significant synchrotron flux from the inner Galaxy. 

An alternative possibility for the origin of the WMAP haze is a single nearby and young 
pulsar, aligned along the line-of-sight toward the Galactic Center. A Vela-like pulsar (r ~ 10^ 
years, d ~ 10^ pc), for example, could yield the observed (radially symmetric) morphology 
and hard spectral index, if located in the desired direction of the sky. To normalize to the 
observed synchrotron flux, however, would require the pulsar to inject ~ 10^^ erg into energetic 
electrons/positrons, making it difficult to understand how such an object would have gone 
undetected in other wavelengths. We thus conclude that this scenario is very unlikely to explain 
the presence of the WMAP Haze. 

V. CONCLUSIONS 

In this paper, we have calculated the intensity and spectrum of synchrotron emission from the 
highly relativistic electrons and positrons produced in relatively young pulsars. In particular, 
we considered the 376 pulsars in the ATNF catalog which are less than 10^ years old, and 
produced maps predicting the resulting synchrotron emission. 

If pulsars inject ~10^^ erg or more into ~10-1000 GeV electrons/positrons, as would be 
required for such an object to generate the rising positron fraction observed by the PAMELA 
experiment, they could also produce bright and potentially observable features in the syn- 
chrotron sky. This is especially true in the case of very nearby and young pulsars (such as the 
Vela or Crab pulsars). In the region of the sky within ~30° of the Vela pulsar, for example, we 
predict considerably brighter and harder synchrotron emission than in other regions of the sky 
away from the Galactic Plane. A smaller region of very bright and hard synchrotron emission 
is also predicted from the region around the Crab pulsar. 

Intriguingly, the features predicted by our calculations are qualitatively similar to those 
observed in the WMAP residual maps (the emission not correlated with other known fore- 
grounds) [32] ■ Given the astrophysically complex nature of these regions (the Gum Nebula and 
Crab Nebula), however, it is not possible at this time to positively identify the origin of this 
emission. 

This work has been supported by the US Department of Energy, including grant DE-FG02- 
95ER40896, and by NASA grant NAG5-10842. 
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